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This work deals with braneworld models driven by real scalar fields with nonstandard dynamics. 
We develop the first-order formalism for models with standard gravity but with the scalar fields 
having generalized dynamics. We illustrate the results with examples of current interest, and we 
find analytical and numerical solutions for warp factors and scalar fields. The results indicate that 
the generalized braneworld scenario is classically stable, and capable of localizing gravity. 
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I. INTRODUCTION 

The appearance of electromagnetism in the nineteenth 
century has triggered a fundamental question, which cul- 
minated with the establishment of Special Relativity. 
The relativity concept is based on the Lorentz group, 
which was born confronting Galilei invariance and its 
wrong concept that information may travel at arbitrarily 
large speed. Not too much later, Special Relativity gave 
birth to General Relativity, and nowadays, almost one 
century later, one is facing another consistency problem, 
this time related to the Planck length. 

Since Relativity makes the speed of light a fundamental 
constant of nature, it necessarily requires that both space 
and time change when the coordinate system is changed 
to another one, with different speed, and this seems to 
confront with the concept of a fundamental length, the 
Planck length. We can add to this problem the apparent 
inconsistency between General Relativity and Quantum 
Mechanics, and the so-called dark energy problem, which 
has appeared due to the recent experimental observation 
that the universe is undergoing accelerated expansion. 

The above issues have currently led researchers to con- 
sider the possibility of including modifications of the 
standard scenario involving matter and geometry. Sev- 
eral possible ways are under consideration, and a very 
popular procedure is known as the quintessence way, in 
which one in general includes dynamical scalar fields that 
can interact through a diversity of possibilities. 

In the present work, we will focus our attention on the 
braneworld scenario, but we will deal with scalar fields 
with nonstandard kinetic terms coupled with standard 
gravity. In the well-known Randall-Sundrum model 
we can further add scalar fields Q with usual dynam- 
ics and allow them to interact with gravity in the stan- 
dard way. In this scenario, the smooth character of the 
solutions generate thick brane with a diversity of struc- 
tures [1, 0, HJ. However, even when one look for static 
solutions, the intrinsic nonlinear character of the Ein- 
stein equations usually result in an intricate system of 
coupled ordinary differential equations that are hard to 
solve. Despite the possible numerical treatment, it is 
also of interest to find models which support analytic 
solutions. In this particular, one can consider specific 



situations where first-order differential equations appear 
describing the scalar field and metric functions, with the 
potential having a very specific form [3|, |4j, |5| . 

In recent years, there appeared some interesting mod- 
els with non canonical dynamics with focus on early time 
inflation or dark energy, as good candidates to solve the 
coincidence problem [6|. These kind of models have also 
been discussed in investigations of topological defects. In 
fact, global topological defects have been considered in 
0> fl S El • For instance, in Ref. Q one has found do- 
main walls, global strings and global monopoles, and in 
Ref. [§[ some formal aspects of unidimensional topolog- 
ical solutions has been studied, and there it was shown 
that the linear stability is preserved for some classes of 
models. Also, in 0] it was shown that quartic poten- 
tial can support compacton solutions for a specific La- 
grangian, and some local vortices were investigated in 
fioj l . In Ref. the generalized models are used for 
stabilization of inter-brane distance. The brane study 
was also considered in two other works [1^ ], under the 
action of specific scalar field model which gives rise to 
compacton solutions. 

The search for analytical solutions for such general- 
ized models is a non trivial task. Our motivation has 
arisen from a previous investigation, in which one con- 
siders first-order differential equations to solve the cor- 
responding equations of motion [13| . The presence of 
first order equations simplifies the investigation and can 
yield to other extensions, as we will show in this work, 
where we modify the standard braneworld scenario with 
the inclusion of scalar fields with nonstandard dynamics. 

We develop the investigations as follows: in the next 
Sec. Ql] we study (4, 1) brane models with nonstandard 
kinetic term coupled with gravity. As usual, we suppose 
that both the scalar field and the warp factor depend 
only on the fifth extra dimension. After an Ansatz for the 
metric characterizing M4 branes with an asymptotically 
AdSc, bulk is implemented, the general structure of the 
equations of motion are obtained. We show from stability 
analysis that metric and scalar field perturbations can 
be decoupled in the transverse traceless gauge and that 
the spectrum of excitations produces a zero-mode and a 
continuum of positive massive modes. This shows that 
the proposed generalized scenario is stable and capable of 
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localizing gravity in a way similar to the standard case. 
Encouraged by this general result, in the next Sec.lIIIIwe 
implement the first-order framework put forward in [l3j , 
in order to investigate two distinct families of models. 
In Sec. IIVI we end the work with some comments and 
conclusions. 

The generalized model that we consider has the form 
C = F(X) — V(4>), and below we study two specific forms 
for the nonstandard kinetic term. The first one is given 
by F = X -\- a\X\X and depends on a parameter a which 
drives the model away from the standard case. Analyt- 
ical expressions for small a are obtained and compared 
with the numerical investigation done for a larger range 
of values of a. We investigate the energy density of the 
brane as well as the necessary conditions for gravity lo- 
calization and modifications of Newton's law, and there 
we show that the numerical study confirms the analytical 
results obtained for small a. Here an interesting result is 
that gravity localization seems to be more effective at 
smaller values of a, showing that the robustness of the 
model seems to weaken for increasing values of a. The 
second model is described by F = —X 2 , and although it 
is well distinct from the former one, our classical investi- 
gation show that it is also capable of localizing gravity. 



II. GENERALIZED DYNAMICS 

We start with a five-dimensional action in which grav- 
ity is coupled to the scalar field in the form 



S 



R + C(6,X' 



(1) 



where we are using = l/(47r) and the signature of 

bulk metric as (H ), with g — det(gMJv)- We 

take the spacetime coordinates and fields as dimension- 
less quantities, and the convention M,N = 0,1,2,3,4 
and fi, v = 0, 1, 2, 3. We also define the invariant 



1, 



X = -V m </>Va/<; 



(2) 



The Einstein equations are Gab — 2 Tab, with the 
energy-momentum tensor having the form 

Tab = V ' a4>V ' b 4> C x - 9abC (3) 

The equation of motion for the scalar field is given by 

G ab V a Vb<P + 2XC X4> -C^ = (4) 



where G AB has the form 

G AB =£ x g AB + £ xx V A <pV B < 



(5) 



We use the notation Cx = dC/dX and Cf, = dC/d(f>, 
etc. In order for the above differential equation to be 
hyperbolic, the condition 



£ 



X 



2XC 



XX 



Cx 



> 



(6) 



must be fulfilled. 

We use the standard notation, and write the metric as 



ds 2 = e 2A X] lLV dx ,1 dx v - dy 2 



(7) 



where A = A{y) describes the warp factor and only de- 
pends on the extra dimension y. As usual, we suppose 
that the field <j) is static, and also, it only depends on the 
extra dimension. Thus, A = A(y) and <f> = 4>(y), and so 
the equation of motion for the scalar field reduces to 

{C x + 2XC X x) <t>" - {2XC X4 , ~ £*) = -4C X <I>'A' (8) 

where prime denotes derivative with respect to the extra 
dimension. The Einstein's equations with the metric (O 
lead to 



A" 
A' 2 



\XC X 

i (C - 2XC X ) 



(9a) 
(9b) 



where for static solutions we have X — —<fi' 2 /2. The 
equations ((SJ and © are not independent, the last one 
being the null energy condition that imposes a brane with 
positive pressure, obeying {C — 2XCx > 0). In particu- 
lar, we can multiply ([8]) by <f>' in order to get 



(C - 2XCx)' = -H' 2 A'C X 



(10) 



and now, if we substitute Eq. (|9b[) we then recover Eq. 
([9lip. 

We note that the Eqs. (|8j)- (|9b[) reduce to the known 
equations in the standard case, in which C = X — V : 

<t>" + Acf) 1 A 1 + Vj, = (11a) 

A" + -<j>'* = Q (lib) 
3 

A 12 - -U' 2 + lv(4>) = (11c) 
6 3 

An important characteristic of the brane is its tension, 
which is given by 



T = J dye 2A ^T 00 = J dy Pl 



(12) 



where p(y) = —e 2AI - y ^C is the energy density. 

The proposed investigation may be of direct interest 
to high energy physics, but it is important to know if 
the modification of the scalar field dynamics will con- 
tribute to destabilize the geometric degrees of freedom of 
the braneworld model. We investigate this issue study- 
ing linear stability in the usual way. We consider metric 
perturbations in the form 



ds 2 



? 2A[V) (V^ + K v (y,x))dx»dx v - dy 2 (13) 



We must also consider fluctuations of the scalar field 



= <t>{v) + ^(y^) 



(14) 
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The first order contribution of the fluctuations to the 
scalar X is written as = (l^h^d^d^ + d^^d^. 
We found the first order contributions of Einstein equa- 
tions in Ricci form as Rab — Tab, with Tab — Tab - 
(l/2>)g AB T% and 

r(!) - -t^p 2 - 4 

-X (Cx<p4> - £xx4>'4>j + 

2e 2A V ( X C X - C) (15a) 

(15b) 



f(i) 



3 ( 2 ^ x< ^ + ^ ^ 
~(2£ XX X + 3£ X )<P'& 



In this case, Einstein's equations turn out to be 
e 2A (l/2d 2 + 2A'd y ) V + ^V« M «'9y(/^) 



-?7 Q/3 {d^duha/3 - dfj.da.Kp - dvdah^p) 



and 



The next steps are known: we introduce the z- 
coordinate in order to make the metric conformally flat, 
with dz = e~ A ( y ^dy and we write 



HUz) = e-**e 3 / 2>1 W V 



(21) 



In this case, the 4-dimensional components of h^v obey 
the Klein-Gordon equation and the metric fluctuations 
of the brane solution lead to Schrodinger-like equation 



where 



[-d 2 + U(z)] = p 2 H, LV 



U(z) = \a*(z) + \a"(z) 



(22) 



(23) 



(15c) We- can write this equation in the form 



d z + \A'{z) \ (-d, + \a'(z) ) = p 2 H,„ (24) 



This factorization directly shows that there are no gravi- 
ton bound-states with negative mass, and the graviton 
zero mode H fll/ (z) cx e^ A ^ is the ground-state of the 
associated quantum mechanical problem. 

This result leads to the important conclusion that 
the modification of the scalar field dynamics does not 
contribute to destabilize the geometric degrees of free- 
dom which appears in the standard braneworld scenario. 
Thus, the modification here proposed is robust and may 
be of direct interest to high energy physics. 



-ifd (dahup - dphap) = £x<$>d^ 
-\(d 2 y +2A' 2 d y )^h aP 

= -^(2£ Xct> X + £^4>+l(2£ xx X + 3£x)4>'4>\l7) 
The equation of motion for the scalar field gives 

£ x e 2A n^-( K (2£ xx X+£ x )^'-AA'{2£ xx X+£ x )& 
- (4£x^'A' + {£ x ^')'+£^) 4> = £x<P'r, aP h' ap (18) 

Let us now consider the transverse traceless compo- 
nents for metric fluctuations 

Kv = ^ (V 1 ^ + n ^f3 7T tJ.a) — ^nn^a/^j h af3 (19) 

where n^v = i]^ iy — d fl d l/ /D. We note that the net effect of 
this projection operation is to decouple the metric fluc- 
tuation equation from the scalar field equation, even in 
the general case which is being considered in the present 
work. Indeed, we can check that equation (|16[) reduces 
to the known equation 



III. THE BRANEWORLD SCENARIO 

Let us start reviewing the case without gravity, setting 
A(y) = 0, which means that only the scalar field equation 
of motion has to be considered. We follow d, [l]| and we 
get 



(Cx + 2X£xx) 4>" = 2X£x<f> - £4 



(25) 



In this case, we can use Eq. (flU)) to get £ — 2£ X X = C, 
where C is an integration constant which can be identi- 
fied with the pressure T in the absence of gravity. For 
stable configurations, the pressureless condition is neces- 
sary. Thus we write 



£ - 2£ X X = 



(26) 



}2 , A A'a e -2A 



(d 2 y +4A'd y 



□) Kv = 



(20) 



This equation depends on the scalar field and its first 
derivative. Therefore, it is a first-order equation. The 
tension of the solution is 



T = - I dy£= I dy£ x cj)' 2 (27) 



If we introduce the function W — W(4>) such that 

Cx<p' = (28) 
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we can write the tension in the form 

T = W(<p(oo)) — W((p(—oo)) (29) 

The interesting thing here is that the tension does not 
depend on the explicit form of the solution, but only on 
its asymptotic values. 

Let us now consider some explicit examples. For in- 
stance, we can deal with 

C = X + a\X\X - V{4>) (30) 

where a is a real, non negative parameter. We name 
this the type I model. Alternatively, we could choose the 
function 

C = X - aX 2 - V{4>) (31) 

but this does not change the classical scenario, as high- 
lighted in Ref.[l3|. Of course, if a = the standard 
scenario is restaured. For a a very small parameter, we 
can investigate the contribution of this term as a pertur- 
bation to the standard scenario. This leads us with the 
expressions 

(/>' = W^-aWl (32) 

V(<f>) = \wj - (33) 

We can also consider another model, for instance 

C = -X 2 - V(<p) (34) 

We name this the type II model, and here the first order 
formalism leads to the equations 

0' = W$ (35) 

V(4>) = jW| (36) 

More details of the first-order formalism in flat spacetime 
can be seen in d, [l3j] ■ 

To extend the first-order framework to the braneworld 
scenario, we follow some works in Ref. and choose the 
derivative of the warp factor with respect to the extra 
dimension to be a function of the scalar field, and we 
write 

A' = ~W(cf>) (37) 

This equation also appear in the standard braneworld 
scenario. Our point here is that since in the action (|TJ) 
the geometric sector remains unchanged, we expect that 
this equation remains unchanged too. In this case, we 
use the equation (|9"a| to write 

C x <t>' = iw* (38) 

which is the same that appears in the absence of gravity. 
The null energy condition (|9b[) can be written as 

C-2C x X = \w 2 (39) 



The equations ((38)) and ([39]) impose a constraint in the 
Lagrange density. It is not difficult to show that these 
first-order equations solve the second-order equation of 
motion ijHJ). 

Let us consider the standard braneworld model. The 
Lagrange density for the scalar field is given by 

C = X - V((j)) (40) 

The set of equations (J3TJ) , {38]) and (39]) give 

4>' = \w* (41a) 
and the constrained potential is 

V{4>) = \wl - \w 2 (41b) 

These equations describe BPS solutions since they appear 
in supergravity [3[. 

Now we consider the case in which the scalar field has 
generalized dynamics. The general structure of the La- 
grange density is given by 

C = F(X) - V(<P) (42) 

For such models the scalar field is sometimes called a 
k-field. The equation of motion is 

(F' + 2XF")(f>" -V 4> = -4F'(j>'A', (43) 

and from the Einstein equations we obtain 

A" = -F'X, A' 2 = \{F —V — 2XF') (44) 

We notice that in the standard situation F(X) — X, the 
above equations lead to the standard braneworld case. 
In the generalized situation, the first-order equations are 
([37]) and 

F'4>' = (45a) 

F — 2F'X — V{4>) = -W 2 (45b) 
3 

The equation (|45aj) has the form G(<p') = \W<f,. For some 
models this can be rewritten as 

0' = G- 1 (46) 

Now, substituting this into Eq. (|45b[) leads to the poten- 
tial 

V{4>) — (F — 2F'X)\^ =a - Hiw ^ - \w 2 (47) 

Comparing this expression with Eq. (|41bj) for the poten- 
tial of the standard case, we see that only the W$ portion 
is changed. The reason for this is that the W portion fol- 
lows from a pure geometric contribution which remained 
unchanged. 

We now illustrate the investigations with the type I and 
type II models, as described in (|3T)]) and ([3~4"|) , respectively. 



A. Type I model 

We use the function F = X + a\X\X . In this case, the 
equation of motion for the scalar field becomes 



0" + A(f>' A' - V* = a (30" - 40'A')0' 2 



(48) 



Using the first-order equations, we see that the scalar 
field has to obey 



1, 



b' + a0' 3 = -Wt 



(49) 



This algebraic equation of third degree for 0' has the only 
real solution 



t(Wg) 2 
6a m(W^) 



(50) 



where 



m(W^) = (54a 2 TU0+6\/3(16a 3 +27a 4 W^|) 1/2 ) 1/3 (51) 
From Eq. (|47|) we can write the potential as 

(52) 

1 o 

or then, substituting ([50)) 

1 /m(W*) 2 



7w 



6a m{W r j > ) 



+ t ( ^ - - 3^ (53) 

In order to ease investigations, let us focus our study 
in the case of a very small. Here we get, up to first-order 
in a, the field equation 



with the corresponding potential 



The equation (|54[) yields, after an integration, 



(54) 



(55) 



(56) 



and so 0(y) = 0o (y — aW(<j)o)/2), where 0o(y) is the 
solution when a vanishes. We expand this solution to 
get 



Hy) = My)-^4>(yW(My)) (57) 



or using Ij53]l 



Hy) = My) - jW^Mv^wiMy)) (58) 



The warp factor obeys the equation 

A' = -(l/3)W(My) - (a/2) #,(yM0o(tf))) (59) 
It is then easy to see that 

A(y)=A (y) + ^W(My)) 2 (60) 

where Aq represents the standard warp factor, for a = 0. 
The brane tension (fT2l is 
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(61) 



or better 



where 7o is the tension of the brane in the standard sce- 
nario. 

We can consider the explicit example, with W(<p) given 
by 



W{4>) = 3asin(60), 



(63) 



i 
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FIG. 1: Plots of the potential J53J for & = v^/3, a = 1 and 
for a = 0.1 (solid line), a = 1 (dashed line) and a = 10 
(dot-dashed line). 

In this case the potential for a small has the form 

Q 81 a 

V = ^a 2 b 2 cos 2 {b<t>) - 3a 2 sin 2 (60) - ^ a 4 6 4 cos 4 (60) (64) 

In Fig. 1, we plot the potential for some values of a, 
not necessarily small. We have also depicted the other 
Figs. 2,3, and 4 below, numerically, for several values of 
a, and we have also checked that the corresponding an- 
alytic expressions, obtained for a very small, completely 
agrees with the numerical results with a small. We then 
note that the numerical study give full support to the 
analytical expressions which we have obtained for a very 
small. We further note that the numerical study shows 
the robustness of the model for a large range of possibil- 
ities for the parameter a. 
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We can use (|58l) to write for the scalar field 4>(y) 
0o (j/) — (9aa 2 6/8) sin(260 o ), or explicitly 



4>{y) 



arcsm 



tanh ( 2 a ^ 2 y 



9aa 2 b , /3 \ , /3 , 
— tanh I - ab y I sech I - ab y ) (65) 



Fig. 2 shows the kink profile (upper panel), with the a 
parameter increasing the brane thickness, as can also be 
seen from the plots of the energy density (lower panel). 
This is an interesting result, since it shows that the 
a-dependent term used to modify the dynamics of the 
scalar field contributes to thicker the brane. 




FIG. 2: Plots of the scalar field 4>{y) (upper panel) and the 
energy density p(y) (lower panel) for a = 0.1 (solid line), 
a — 1 (dashed line) and a = 10 (dot-dashed line). We use 
b 2 = 2/3. 

We use Eq. dSDJ) to get 

A(y) = ^ln(scch g a 6 2 y ))+^tanh 2 Qa6 2 y )(66) 

which is depicted in Fig. 3. For b 2 — 2/3, we use the 
transformation dz = e~ A ^dy to write 

3a 2 

dz — cosh(ay) — sinh (ay)sech(ay) (67) 

where 

sinh(ay) 3aa , . . , . . , . 

z = 1 — — (arctan(sinh(ay)) — smh(ay)) (68) 




FIG. 3: Plots of the warp factor e 2A ^ (upper panel) and 
the Schrodinger-like potential U(z) (lower panel) for a — 0.1 
(solid line), a = 1 (dashed line) and a = 10 (dot-dashed line). 
We use b 2 = 2/3. 



y = — arcsinh(az) 
a 



3aa (arctan(az) — az) 



(69) 



4 VI + a 2 z 2 
The warp factor is now written in terms of the z variable 



A(z) 



1 



ln(l + a 2 z 2 ) + 



3a 3 a z arctan(a2:) 



2 ~ x ' " ' ' 4 1 + a 2 z 2 
The Schrodinger-like potential is given by 

U{z) = \a*(z) + \a»{z) 

and has the explicit form 

3a 2 (5o 2 z 2 -2) 9 a 4 a 



(70) 



(71) 



U{z) 



4(l + a 2 z 2 ) 2 8(i + a 2 z 2)3 
(az(5a 2 z 2 - 9) arctan(az) + (2 -7a 2 z 2 ))(72) 



In Fig. 3 (lower panel) we plot the potential U(z) for 
several values of a. We see that the characteristic volcano 
profile for dynamically generated M4 brane immersed in 
an asymptotically AdS$ bulk. The increase of a leads to 
a reduction of the maximum of the potential, and this 
may modify the way gravity is localized in the brane. 
Thus, we have to investigate if this behavior produces 
any sensible effects on the localization of gravity in the 
brane generated in this case. 



7 




5 



FIG. 4: Plots of z 2 U(z) for the type I model (upper panel) 
for a = 0.1 (solid line), a = 1 (dashed line) and a = 10 
(dot-dashed line), showing the asymptotic regime for large 
z. We use b 2 = 2/3. We also plot the ratio Ml /Ml (lower 
panel) between 4-dimensional coupling and fundamental 5- 
dimensional Planck scale as a function of a. 



To investigate this possibility, let us look at the asymp- 
totic behavior of U(z). In Fig. 4 (upper panel) we plot 
z 2 U(z) 1 indicating that we have U(z) ~ 1/z 2 for large 
z, even for large values of the parameter a. We see 
from this figure that the asymptotic regime is better 
achieved for smaller values of ev, whereas for larger val- 
ues one must consider still larger values of z. For in- 
stance, for —30 < z < 30, Fig. 4 shows that we can 
study analytically the correction for the Newtonian po- 
tential in the range < a < 0.1. This is the region 
where our analytic expansion for small a agrees sensi- 
bly with the numerical simulations. We find U(z) ~ 
(15/4)/z 2 + (45/16)7raa/z 3 + 0(1/ z 5 ) for large z. Note 
that the leading term of the expansion does not depend 
on a. It is well known that potentials which asymptotes 
as U(z) ~ (3((3 + l)/z 2 gives a correction for the New- 
tonian potential 0(1/ R 2 ") for two massive objects at a 
distance R from each other. In our case, we have (3 = 3/2, 
and this gives the correction C(l/i? 3 ) for the Newtonian 
potential, independent of a. This is the same correction 
given by at the standard Randall-Sundrum scenario, and 
it confirms that the model localizes gravity for all chosen 
values of a where the expansion applies. The simula- 
tions also indicate that the same applies for larger values 



of a. This means that a whole class of braneworld models 
were constructed with different properties related to mat- 
ter distribution and geometry, all being able to localize 
gravity. 

Another point is that the next to leading term depends 
on a for the Schrodinger potential. It shows that for 
larger a the asymptotic regime is achieved for larger val- 
ues of z, as demonstrated by the simulations shown in 
Fig. 4. In order to better see the influence of a in the 
gravitational interaction we remind that the Newtonian 
potential are corrected by the contribution of all massive 
modes, solutions of the Schrodinger-like equation, as 

i i r°° p ~mB, 
U{R) = G R + IPj dm ^\^m(0)\ 2 , (73) 

where the 4-dimensional coupling is G = M^ 2 , and M* 
is the fundamental 5-dimensional Planck scale, and the 
integration is considered at the brane position z = 0. We 
can write an expression relating the two scales as 

/+oo 
dze 3A ^. (74) 
-oo 

In this way one can see that those scales are related to the 
integral of a function depending on the warp factor. In 
Fig. 4 (lower panel) we plot the ratio Mf /M 3 between the 
two scales. Note the greater importance of M4 for larger 
values of a. Since G = 1/M 2 , this means that smaller 
values of a contribute to increase the intensity of the 
gravitational interaction, and that gravity localization is 
then favored. From Eq. (|70p we can write M|/M 3 = 
2/a + aa 2 /3 for small a. This is related to the fact that 
the Schrodinger potential achieves the asymptotic region 
at smaller values of z, for smaller values of a. This effect 
seems to be in perfect agreement with the former result, 
which has shown that the brane thickness decreases with 
decreasing values of a, since in a thicker brane, gravity 
localizes less importantly. 

B. Type II model 

We now study models with the function F — —X 2 /2. 
In this case, the equation of motion of the scalar field is 

j^'V + 20' 3 A' - = (75) 

Using the first order formalism, the equation for <j>' is 

4>' = W* (76) 
and the potential has the form 

V(4>) = \W$ - \w 2 (77) 
We choose an explicit function 

W{4>) = 9a 3 b 2 sin(60)(2 + cos 2 ^)) (78) 
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Here the potential has the explicit form 

V(<P) = =^a 4 b 4 cos 4 (&</>) 
8 

-27a 6 b 4 sm 2 (b(f>) (2 + cos 2 (60)) 2 (79) 

In Fig. 5 we plot the potential for a = 1 and several val- 
ues of 6, and for 6=1 and several values of a. We note 
that in the first case, for a = 1 the effect of increasing b is 
to deepen and narrow the potential wells. In the second 
case, however, for b = 1 and for a increasing, each poten- 
tial well deepens and widens, and this nicely contributes 
to distinguish the two cases, as we show below. 




FIG. 5: Plots of the potential for a = 1 (upper panel) and 
for b — 1 (solid line), b — 1.1 (dashed line) and b = 1.2 (dot- 
dashed line), and for b = 1 (lower panel), and for a = 1 (solid 
line), a = 1.1 (dashed line) and a = 1.2 (dot-dashed line). 

The scalar field obeys 4>' = 3ab cos(6</>) , with the solu- 
tion 



4>(y) = — arcsin(tanh(3a6 2 ?/)) 



(80) 



We can also get for A 

1 2 
A(y) = --a 2 tanh 2 (3ab 2 y) + -a 2 In (sech (3 a6 2 y))(81) 

D o 

which results in the warp factor depicted in Fig. 6, for 
b = 1, with a = 1, 1.1, and 1.2; we are not showing the 
other plots, for a = 1 and for b varying, because they are 
essentially the same. 




FIG. 6: Plots of the warp factor for b — 1 and for a = 1 (solid 
line), a = 1.1 (dashed line) and a — 1.2 (dot-dashed line). 
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FIG. 7: Plots of the Schrodinger-like potential f/(z) for type 
II model. In the upper panel a — 1 and b = 1 (solid line), 
b = 1.1 (dashed line) and b = 1.2 (dot-dashed line). In the 
lower panel b = 1 and a = 1 (solid line), a = 1.1 (dashed line) 
and a = 1.2 (dot-dashed line). 



In the plots of Fig. 7, we see that gravity localization 
seems to be favored in the case 6=1, for increasing a, 
since there we see that the height of the maxima of U(z) 
is higher then in the other case, with a — 1 for increas- 
ing b. However, to check this behavior quantitatively, 
we repeated the procedure included in Sec. IIII Al in order 
to numerically obtain the functions z(r),A(z) and U(z). 
The function U{z) was graphically analyzed to investi- 
gate the influence of the parameters a and b for gravity 
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localization. This is done in Fig. [3 in which one displays 
U(z) for several values of parameters a and 6, and there 
we note that: i) both figures have similar volcano profile. 
However, the numerical analysis shows that fixing 6 = 1 
with 1 < a < 3 leads to larger extremes for the poten- 
tial than fixing a = 1 with 1 < b < 3. In this way, we 
can infer that the parameter a has greater influence for 
gravity localization; ii) fixing one of the parameters a or 
6, the increasing of the second parameter leads to a cor- 
responding increasing on the maxima and minima of the 
potential, indicating that the increasing the parameters 
favors gravity localization. This is better seen in Fig. [51 
where we study z 2 U(z). 




FIG. 8: Plots of z 2 U(z) for the type II model. In the upper 
panel a — 1 and b = 1 (solid line), 6 = 2 (dashed line) and 
b = 3 (dot-dashed line). In the lower panel 6=1 and a = 1 
(solid line), a = 2 (dashed line) and a = 3 (dot-dashed line). 

In Fig. [8jone notes that the behavior of U(z) suggests 
that U(z) ~ 1/z 2 for z » 1. Writing U(z) = /3(/3+l)/z 2 
we can, as in the analysis of the type I model, determine 
(3 such that the Newtonian potential has correction of 
order l/i? 2/3 for large separation R between unit masses. 
Now from these figures we can conclude that: i) for a = 1 
and fixed z max , the asymptotic approximation is easier 
obtained for larger values of b. On the other hand, larger 
values of b leads to larger values of f3, and this corre- 
sponds to a correction with a greater power law for the 
Newtonian potential. This agrees with our previous anal- 
ysis that larger values of a or b favor gravity localization; 
ii) for z max — 20 and a = 1, we can affirm that the 
asymptotic region where U(z) — ({3)((3+l) / z 2 is achieved 



with considerable precision for b > 2. From the numer- 
ical values obtained we can estimate that: for 6 = 2, 
we get (/3)(/3 + 1) = 3.7180, (3 = 1.492, and a correc- 
tion C(l/i? 2 - 98 ) for the Newtonian potential (compare 
with 0(l/i? 3 ) for Randall-Sundrum model); for 6 = 3, 
we get (/?)(/? + 1) = 3.7338, (3 = 1.496, and a correction 
0(1/ R 2 ") for the Newtonian potential. 

We further note that for z max = 20 and 6 = 1 we can 
affirm that the asymptotic region is achieved for a > 2. 
With the same procedure we estimate that: for a = 2, 
we get (/?)(/? + 1) = 3.7809, = 1.508, and a correction 
C(l/i? 3 02 ) for the Newtonian potential; for a = 3 we 
get ((3)(f3 + 1) = 3.8565, (3 — 1.526, and a correction 
0(l/R 3 05 ) for the Newtonian potential. 

The above comments agree with the previous obser- 
vation that larger values for a and 6 lead to corrections 
for the Newtonian potential with larger power law. The 
corrections are near to C(l/i? 3 ) obtained from Randall- 
Sundrum model for small parameters, and tend to be 
larger for larger a and 6, showing that gravity is then 
easier localized. 

If we compare the case a — 1,6 = 3 (correction 
0(l/R 2 ")) with the case a = 3,6 = 1 (correction 
C(l/i? 3 05 )), we see that the increasing of the parameter 
a is more effective for gravity localization in comparison 
with similar increasing of parameter 6. 



IV. ENDING COMMENTS 

In this work we have investigated generalized 
braneworld models, that is, models in which the stan- 
dard gravity action is extended to include scalar fields 
with generalized dynamics, with the Lagrange density 
having the nonstandard form C(4>,X) = F(X) — V ((/>). 
This study is a continuation of our program to investigate 
the scalar field behavior under the presence of generalized 
dynamics [l,[l3j]. In particular, in the present study we 
have included two distinct families of models, one given 
by F(X) = X- aX\X\, and the other by F(X) = -X 2 . 

An interesting and general result is that the proposed 
scenario, in which gravity acts standardly and the scalar 
field has generalized dynamics, is linearly stable, capable 
of localizing gravity in a way similar to the standard case. 

Other investigations included in this work engender 
both analytical and numerical results. For the type I 
model, with F(X) = X — aX\X\, we have presented an- 
alytical results up to first order in a for a very small, and 
our numerical study confirms the correctness of the per- 
turbative expansion up to first order in a. In this case, 
for a small, we have calculated analytically the correc- 
tions to the Newtonian potential, showing the localiza- 
tion of gravity. The ratio between the masses leads to a 
nice estimate of how a quantitatively affects the gravi- 
tational interaction. And numerically, we could extend 
this result to much larger values of a. As another inter- 
esting result, we have shown that gravity localization is 
more effective at smaller values of a. It seems that for 
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F(X) = X + a\X\X , the robustness of the model weak- 
ens for increasing values of a, as we get away from the 
standard braneworld scenario. 

The type II model is more involved. However, thanks 
to the first-order framework put forward in [131 ] , we could 
study it and obtain analytic solutions for both the scalar 
field and warp factor. The analytic solutions has helped 
us to ease the subsequent numerical investigations, to 
study gravity localization. The results show that, like 
in the former case, this new and well distinct family of 
models engenders similar behavior, and it is also capa- 
ble of localizing gravity. For the two families of models, 
we can also control gravity localization with the specific 
form of the potential V(<fi), an effect that also appears in 
the braneworld model with the scalar field with standard 



dynamics. 

There are other possibilities of study. For instance, in 
the case of the type II model, we can find compacton 
solutions for the scalar field, so we can also investigate 
this case, in a way similar to the study done in [121 ]. 
Another possibility is to extend the present investigations 
to the case of bent brane, as we have already done in some 
of our work in Q in the case of standard dynamics. We 
can also consider the harder case, in which we also change 
R to F(R), extending the F(R) brane study done in [14[ 
to this new F(R, X) braneworld scenario. We hope to 
report on these issues in another work in the near future. 

The authors would like to thank CAPES, CNPq, 
CNPq-MCT-CT-Encrg and PRONEX-CNPq-FAPESQ 
for partial support. 
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